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Abstract. We study singularities of the n-body problem in spaces 
of constant curvature and generalize certain results due to Painleve, 
Weierstrass, and Sundman. For positive curvature, some of our 
proofs use the correspondence between total collision solutions of 
the original system and their orthogonal projection — a property 
that offers a new method of approaching the problem in this par- 
ticular case. 

1. Introduction 

We consider the n-body problem in spaces of constant curvature, 
which we will henceforth call the curved ra-body problem to distin- 
guish it from its classical Euclidean analogue. Our goal is to study 
solutions that experience total collisions and, to some extent, solutions 
that end in some kind of hybrid singularities, i.e. both coUisional and 
non-collisional. 



1.1. Results. We generalize in Section 3 a criterion proved by Paul 
Painleve in 1897, which shows that a solution of the Euclidean n-body 
problem has a singularity if and only if the limit of the minimum dis- 
tance between particles tends to zero, [12]. Our generalization takes 
into account some singularities with no correspondent in the Euclidean 
case. But if we disregard them, a step we show to be natural from the 
physical point of view, our generalization reflects Painleve's original 
result. This property, as well as other results obtained in our previous 
work [4], support the more than a century old idea that the potential 

1 



2 



Florin Diacu 



given by the cotangent of the distance is the correct generahzation of 
the Newtonian gravitational model to spaces of constant curvature. 

In the second part of this paper, we establish a connection between 
total collision solutions and the integrals of the angular momentum in 
the case of positive curvature. Two classical results of the Euclidean 
case, the first already known to Karl Weierstrass in the 1880s, were 
proved by Karl Frithiof Sundman in an articl^ published in 1912: (i) 
if a solution of the n-body problem experiences a simultaneous total 
collision, all three angular momentum constants are zero, and (ii) if 
a triple collision occurs in the 3-body problem, the motion is planar. 
We will show in Sections 5 and 6 that, in a suitable setting and under 
some reasonable assumptions, these properties have analogues in spaces 
of positive constant curvature. The proofs use a lemma (developed in 
Section!!]), which establishes, under certain restrictions, the equivalence 
between the equations of motion and their orthogonal projection in the 
Euclidean space of the same dimension as the original phase space. 

1.2. Some historical remarks. The study of gravitation outside the 
Euclidean context started in the 1830s, when Janos Bolyai and Nikolai 
Lobachevsky extended the Newtonian 2-body problem to the hyper- 
bolic space. These co- discoverers of hyperbolic geometry independently 
proposed a gravitational force proportional to the inverse area of a 2- 
dimensional sphere having the same radius as the distance between 
bodies. Ernest Schering showed in 1870 that the terms of the poten- 
tial involve the hyperbolic cotangent of the hyperbolic distance, [T3] . 
Some years later, Wilhelm Killing adopted the cotangent potential in 
the positive curvature case too, [B]. Subsequent studies of the 2-body 
problem proved that Kepler's laws admit natural generalizations, [7|, 
as does Bertrand's theorem, according to which there exist only two 
analytic central potentials for which all bounded orbits are closed, [8]. 
Other attempts (such as those of Rudolph Lipschitz, [9]) at generaliz- 
ing the problem using potentials that do not involve the cotangent of 
the distance failed to recover the classical properties of the Euclidean 
case. A more detailed history of these developments is given in [1] , the 
first piece of work that derives and studies the equations of motion for 
any number of bodies, and proves Saari's conjectur^ in the geodesic 
case. Paper [1] is also a prerequisite for appreciating the results we 
prove in this article. 

"'^This famous 1912 paper, |17] . was in fact an invited exposition to Acta Mathe- 
matica of two previous research articles Sundman had published in 1907, |15] . and 
1909, |16| . in an obscure Finish journal called Acta Societatis Scientiarum Fennicae. 

^Details about Saari's conjecture can be found in [5] and [3]. 
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2. Equations of motion 

We first describe tlie equations of motion on 2- dimensional manifolds 
of constant curvature, namely spheres embedded in for k, > and 
hyperboloid^ embedded in the Minkovski space M'^ for k < 0, and will 
discuss in Section the generalization to three dimensions for k > 0. 

Consider the masses mi , . . . , m„ > in M'^ for k, > and for 
K < 0, whose positions are given by the vectors = {xi,yi,Zi), i = 
1, . . . , n, and let q = (qi, . . . , q„) be the configuration of the system. 
We define the gradient operator with respect to the vector qj as 

where 



cr 



+1, for K > 
-1, for K < 0, 



and let V denote the operator (Vq^, . . . , Vq„). For the 3-dimensional 

vectors a = {ax, ay,az) and b = {bx,by,bz), we define the inner product 

(1) a © b := {a^bx + ayby + aa^b^) 
and the cross product 

(2) a (g) b := {ttyb^ - a^by, aj)^ - aj)^, (y{aj)y - ayb^)). 

The Hamiltonian function of the system describing the motion of the 
n-body problem in spaces of constant curvature is 

i/«(q,p)=T.(q,p)-t/.(q), 

where 



n 



is the kinetic energy and 



i=l 



"'The hyperboloid corresponds to Weierstrass's model of hyperbolic geometry 
(see Appendix in [i]). 
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is the force function, —11^ representing the potential energjl^. Then the 
Hamiltonian form of the equations of motion is given by the system 



(4) 



Vq,f/K(q) - ^K(pi Pi)qi, i = l,...,n, k 0, 



the gradient of the force function having the expression 
(5) 

V ri ( ) = mimj{aKy/^{Kq^j Q qj)[(/tq.^ q»)qj- - {ncji Q qj)qi\ 

The motion of the bodies is confined to the surface of nonzero con- 
stant curvature k, i.e. (q, p) G T*(M^)", where 

= {(x, y, z)eM.^ \ k{x^ + y^ + (yz^) = 1} 

(in particular, = and M^^ = H^) and T*(M^)" is the cotangent 
bundle of the configuration space (M^)". For k > we will also denote 
by S^, while for k < we will denote it by H^. 
Notice that the n constraints given by Kq^ © = 1 imply that 
qi © Pi = 0, so the 6?7,-dimensional system (jlj) has 2n constraints. The 
Hamiltonian function provides the integral of energy, 

^K(q,p) = h, 

where h is the energy constant. Equations (j4]) also have the integrals 
of the angular momentum. 



(6) ^q.i^p 



i=l 



where c = 7) is a constant vector. Unlike in the Euclidean case, 
there are no integrals of the center of mass and linear momentum. Their 
absence complicates the study of the equations because many of the 
standard methods don't apply anymore. The curved ra-body problem 
is thus a fresh source for new mathematical developments. 

3. Singularities 
Equations (jlj) are undefined in the set A := Ui<j<j<„Ajj, with 
A,, := {qe (M^r I («:q,©q,.)2 = l}. 



^In [3] , we showed how this expression of U^, follows from the cotangent potential 
for K 0, and that Uq is the Newtonian potential of the Euclidean problem, 
obtained as k ^ 0. 
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where the force function and its gradient have zero denominators. Thus 
A contains all the singularities of the equations of motion. The singu- 
larity condition, (Kq, q^)^ = 1, suggests that we consider two cases, 
so we write Aj^ = A^^ U A^:^., where A^^ := {q G (M^)" | /tq^ ©q^ = 1} 
and A~. := {q G (M^)" | nqi Q qj = —1}. Accordingly, we define 

A+ := Ui<i<j<nA+. and A" := Ui<i<j<„A-.. 

Then A = A+ U A~. The elements of A"*" correspond to collisions 
for any k, ^ 0, whereas the elements of A~ correspond to antipodal 
singularities for n > 0. The latter occur when two bodies are at the 
opposite ends of the same diameter of a sphere. For n < 0, antipodal 
singularities do not exist. 

The set A is related to singularities which arise from the question 
of existence and uniqueness of initial value problems. For initial con- 
ditions (q,p)(0) e T*(M2)" with q(0) ^ A, standard results of the 
theory of differential equations ensure local existence and uniqueness 
of an analytic solution (q, p) defined on some interval [0, t~^). Since the 
surfaces are connected, this solution can be analytically extended 
to an interval [0,t*), with < t+ < t* < cx). If t* = oo, the solution is 
globally defined. But if t* < cxd, the solution is called singular, and we 
say that it has a singularity at time t*. 

We have seen in [1] that, for k > 0, no solutions encounter antipodal 
singularities alone. But there are solutions that encounter collision sin- 
gularities and solutions that encounter collision-antipodal singularities, 
as for instance when two bodies collide at the north pole while a third 
body is at the south pole. 

3.1. Generalization of Painleve's theorem. A classical result due 
to Paul Painleve, in the Euclidean case, shows that an analytic solution 
defined on [0,t*) has a singularity at t* if and only if the inferior limit 
of the minimum of the mutual distances vanishes when t ^ t*, [T2] . 
a detailed presentation of which can be found in pQ. Under a certain 
assumption, we can translate this property to spaces of constant curva- 
ture. To prove these results, we start with a couple of lemmas, which 
generalize some properties known in the traditional literature, [18]. 

Lemma 1. // (q, p) is an analytic solution of equations (jlj), defined 
on [0,t*), with t* a singularity, then 

liminf min K^qj © q,,)^ — 1| =0. 

t^t* l<i<j<n ■' 
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Proof. Notice first that using the constraints Kq^Oqi = 1, i = 1, . . . , n, 
the integral of energy becomes 



" " / M /2 ^ 

■mimj[aK) ' Kcii q 



(7) (P.0P.)-E E 



, 2h. 



Also, from the equations of motion (jlj) we can conclude that 

(ffK)3/2[q. - (Kq^Q< 



mim^- ((Tk)3/2 [q. - (Kqi©q )qi] _i , . 



We can now prove the necessity of the condition in Lemma [TJ For this, 
assume that there is a constant c > such that 

liminf min K^q^ q,)^ — 1| > c. 

t— >t* l<i<j<n 

Then there is a time to in [0, t*) and constants bij > c, with 1 < i < j < 
n, for which [cr — cr(Kqj0q,)^]^/^ > bij for all t in [to,t*). Consequently 
equation ([7]) implies that Yl^=i "^i~^(Pi © Pi) is bounded, so every term 
of this sum is bounded as well. Then equations ([H]) guarantee the 
existence of n constants Si > such that |pj| < Si, i = 1, . . . ,n, and 
therefore all q^ are bounded. 

Writing the configuration vector q as a Taylor series about to with 
integral remainder, 

q(t) = q(to) + {t- to)q(to) + [ {t - r)q(r)rfr, 

J to 

and using the fact that q is bounded, we can conclude that there is a 
vector (q*,p*) in phase space such that limt_»t. (q, p)(t) = (q*,p*), i.e. 
the position vectors q^ and momentum vectors pi have limiting posi- 
tions q* and p*, respectively. So [a — a{Kq* Q q*)^]"*^/^ > bij, therefore 
K^q* © q*)^ — 1| > 0, which means that the distance between particles 
is not zero, and they are neither at a collision-antipodal singularity if 
K, > 0. But then the particles can keep moving. This physical conclu- 
sion suggests a contradiction with the hypothesis by showing that the 
solution is analytic at t*. To prove this fact rigorously, notice that the 
domain of the solution (q, p) depends on the constants Si, therefore on 
bij, and finally on c, but is independent of the initial conditions. So 
by choosing the initial data, along the same solution, close enough to 
t* and (q*,p*), the solution remains analytic at t*, a conclusion which 
contradicts the existence of the constant c as described above, and thus 
completes the proof. □ 
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Lemma 2. Assume that (q, p) is an analytic solution of equations 
dl]), defined on [0,t*), that is bounded away from collision- antipodal 
configurations if k,> 0. Then, if 

liminf min iKOj q,- — 1| = 0, 

t^t* l<i<j<n 

t* is a singularity of the solution. 

Proof. Assume that liminft^t. mini<j<j<„ I Kqj0qj — l| =0. Obviously, 
if q becomes unbounded as t — > t*, then t* is a singularity, and Lemma 
[2] is proved. We consequently assume q to be bounded. Then p is 
bounded as t — t*, so the momentum p is bounded as well. Therefore 
we can conclude from equation ([7j) that U^iq) is bounded as t — t*. 

Recall, however, that the terms defining 11^ have denominators of 
the form [a — a^nq^i © q^YY^"^, so when the quantities a — a^nqi Qcy)^ 
are small, the corresponding terms of the force function become large 
in absolute value. But these terms have one sign if Kqj © q^ is near 1 
and the opposite sign if it is near —1. As we excluded from our hy- 
pothesis solutions that come close to collision- antipodal configurations, 
the quantities Kqj © q^ are bounded away from — L Consequently, as 
liminff^f. mini<i<j<„ l^qi © q,- - 1| =0, limsupi^^^, f/^(q(t)) = oo, 
a contradiction with the conclusion drawn at the end of the previous 
paragraph. So the condition assumed in Lemma [2] makes t* a singular- 
ity. This completes the proof. □ 

Remark 1. The reason for having to exclude collision- antipodal con- 
figurations from the hypothesis of Lemma [2] is connected to a property 
proved in Theorem 1 (iii) of [4]. We showed there that there are choices 
of masses and initial conditions for which a 3-body problem taking 
place in can have finite forces and velocities at a collision-antipodal 
configuration; in other words the solution remains analytic at t*. For 
instance, this is the case when two bodies of mass 4m and a third body 
of mass m move on a great circle of S^, forming at each moment an 
isosceles triangle, and such that the larger bodies collide in finite time 
while the smaller body reaches the diametrically opposed side of the 
circle. So there are orbits that do not experience a singularity at t* but 
for which liminf f^f* mini<j<j<„ l^qj © q,- — 1| = 0. 

We can now state and prove a generalization of Painleve's theorem to 
spaces of constant curvature. Like for the above two lemmas, we split 
the result into two statements, one proving necessity and the other 
sufficiency. 
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Theorem 1. // (q, p) is an analytic solution of equations (jlj) defined 
on [0,t*), with t* a singularity, then 

lim min |(fi;qj q,)^ — 1| =0. 

t^t* l<i<j<n ■' 

Proof. To prove the necessity of the condition, assume that there is a 
constant c > such that 

hmsup min |(/tqi q^)^ — 1| > c. 

Then there exists a sequence of times, (t„)„gN (where N represents the 
set of positive integers), with t^ t* such that | [Kqj(t„)0qj(t„)]^ — 1| > 
c > for all i,j with 1 < i < j < n. This means that there is a positive 
constant b for which t/(q(t„)) < b for all n G N. Equation ([7]) leads 
therefore to the conclusion that Yl^=i ^T^i.Pii'^n) Pj(^n)) < 2(6 + h) 
for all n e N. Therefore there is a constant a > for which |p(tn)| < « 
for all n G M. But as we already showed in the proof of Lemma [H the 
domain of the solution is independent of the choice of initial conditions. 
Consequently we can choose some initial data t„Q, q(tno), p(^rio) with tno 
close enough to t* to make the solution analytic at t*, which means that 
t* is not a singularity. This contradiction proves the necessity of the 
condition, and completes the proof. □ 

Theorem 2. Assume that (q, p) is an analytic solution of equations 
(jlj), defined on [0,t*), that is bounded away from collision- antipodal 
configurations if k,> 0. Then, if 

lim min l^qj q,- — 1| =0, 

t-^t* l<i<j<n 

t* is a singularity of the solution. 

Proof. This result is an obvious consequence of Lemma [2J □ 

Remark 2. As long as, for k > 0, solutions stay away from collision- 
antipodal singularities, the condition in Theorem 1 can be reduced to 

lim min l^qj q,- — 1| = 0, 

t^t* l<i<j<n ^ 

and Theorems 1 and 2 become accurate translations of Painleve's orig- 
inal result to spaces of constant curvature. For k < 0, the accurate 
translation is satisfied without restrictions. 

3.2. Remarks on the nature of singularities. In the classical n- 
body problem solutions can experience two kinds of singularities: col- 
lisions and pseudocollisions. The latter, which occur when the motion 
becomes unbounded in finite time, have been conjectured by Painleve 
for n > 3, [12], and proved to exist for four or more bodies, [10], [19], 
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[5]. We do not know whether such singularities show up in the curved 
n-body problem. The compactness of seems to exclude them for 
K > 0, but they may exist for k < 0. The difficulty of solving this 
problem is compounded by the lack of integrals of the center of mass, 
which played a crucial role in proving the existence of pseudocoUisions 
in the classical case. 

The collision-antipodal singularities of the curved problem raise an- 
other question. Are they due to the coordinates or the potential? The 
physical remarks below, point at the potential. Moreover, shifting the 
center of the coordinate system away from the center of the sphere does 
not remove these singularities. But, of course, this doesn't exclude the 
possibility of finding singularity- free coordinates in the future. 

3.3. Some physical remarks. The antipodal and the collision-anti- 
podal singularities seem to obstruct the natural translation of the dy- 
namical properties of the n-body problem from k = to k > 0. To 
better understand this issue, let us first compare how the force function 
and its gradient vary in the Euclidean and in the curved case. 

Let us start with the 2-body problem. The Euclidean force function, 
^o(ci) = — c[2|; is infinite at collision and tends to zero when 

the distance between bodies tends to infinity. The norm of the gradient, 
I V?7o(q)|, has a similar behavior, which agrees with our perception that 
the gravitational force decreases when the distance increases. 

The behavior of the curved force function ([3]) and the norm of its 
gradient ([5]), however, depend on the sign of k. For n < 0, things are 
as in the Euclidean case. For k > let's assume that one body is fixed 
at the north pole. Then ^7^, ranges from -|-oo at collision to —00 at the 
antipodal configuration, with when the second body is on the equator. 
The norm of the gradient is -|-oo at collision, and the smaller the farther 
the second body lies from collision in the northern hemisphere; it takes 
a positive minimum value on the equator, and is the larger the farther 
the second body stays from the north pole while lying in the southern 
hemisphere; finally, the norm of the gradient becomes +00 when the 
two bodies are at antipodes. 

This behavior of the gradient seems to agree with our understanding 
of gravitation only when the second body doesn't leave the northern 
hemisphere, but not after it passes the equator. In a spherical universe 
with billions of objects ejected from a Big-Bang that took place at the 
north pole, all the bodies would now still be in the northern hemi- 
sphere. But when the boundary of the expanding system approaches 
the equator, many bodies come close to antipodal singularities, so the 
potential energy becomes positive, thus having the same sign as the 
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kinetic energy. Then, by the integral of energy, the potential energy 
cannot grow beyond the value of the energy constant, which, when 
reached, makes the kinetic energy zero. Consequently the system stops 
moving. The motion then reverses from expansion to contraction, in 
agreement with the cosmological scenario of general relativity. So in 
a highly populated spherical universe, the motion is contained in the 
northern hemisphere, away from the equator. Therefore we can restrict 
the study of the case k > to a hemisphere without equator and ignore 
solutions that reach collision-antipodal singularities. 



The above physical remarks show that while the natural setting of 
the case k < is the entire upper sheet of the hyperboloid, we can 
restrict the study of the case k > to the northern hemisphere. The 
rest of this paper is about the positive-curvature case. 

We further introduce the equations of the orthogonal projectioiJl, 
which we will henceforth call the orthogonal system, and will show 
that, under certain circumstances, there is a one-to-one correspondence 
between its total-collision solutions and those of the original equations 
of motion. The orthogonal system, which will help us understand cer- 
tain properties of the original equations, has the advantage of being 
defined in an Euclidean disk, though not endowed with the standard 
Euclidean distance. 

Let us also mention that the Principal Axis Theorem (see [1]) allows 
us to use the orthogonal transformations of the sphere to keep, in some 
suitable basis, the original form of the equations of motion. There- 
fore, without loss of generality, we can always apply the orthogonal 
projection on the xy plane. Similarly, we introduce no restrictions by 
assuming that the total collision we study in the next sections takes 
place at the north pole, (0, 0, k"^''^), of S^. 

Let qj = {xi.Ui) be the orthogonal projection of = {xi,yi,Zi) 
with the constraint Kqj ■ qj = 1, i = 1, . . . ,n, on the xy plane. The 
momenta of the projected variables are pi = mjqi, and the problem 
is restricted to the disk Kq^ ■ q^ < 1. Then we obtain the orthogonal 
system associated to equations (jlj) by dropping the Zi variables from 
these equations. In other words, the orthogonal system has the form 



4. The orthogonal system 



(9) 




^In cartography, this projection is called orthographic, and was already men- 
tioned by Hipparchus in the 2nd century B.C. 
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The orthogonal transformation, however, introduces new (artificial) 
singularities in system They occur when two or more bodies reach 
the same diameter of the disk, a position for which at least one de- 
nominator in Vq-f/^lq) vanishes. We call this a diameter singularity of 
system (Q. In terms of (jlj), diameter singularities correspond to non- 
singular configurations for which two or more bodies reach one of the 
geodesies that pass through the north pole (0, 0, k"^^"^). We call them 
pole-geodesic configurations of system (jlj). Since we are concerned with 
total-collision orbits, these artificial singularities force us to exclude 
from our treatment only a negligible set, as we will further see. 

Consider now the set C of solutions of system (jlj) that encounter a 
total collision at time t* at the north pole, and are free of pole-geodesic 
configurations in some interval [tQ,t*), where to is solution dependent. 
Then C contains all total-collision orbits of system (jlj), except for a 
lower-dimensional set. Indeed, each solution is analytic in [to,t*). So 
should the set of pole-geodesic configurations of a solution of system (jlj) 
have an accumulation point, then by the identity theorem of analytic 
functions the entire solution must be pole-geodesic, i.e. two or more 
bodies move all the time such that, at every instant, they lie on the 
same geodesic passing through the north pole. These constraints define 
the lower- dimensional set we must exclude from C. 

Similarly, we consider the set C of solutions of system (j9j) that en- 
counter a total collision at time t* at the origin of the disk, and are 
defined on some interval [to,t*), with to solution dependent. In other 
words, the orbits of C do not encounter diameter singularities in [to, t*). 

We can now state the following result, which proves the equivalence 
between the sets C and C. 

Equivalence Lemma. Consider the set C of solutions (q, p) of system 
(jlj) that encounter a total collision at time t* and are free of pole- 
geodesic configurations in some interval [to,t*). Also consider the set 
C of solutions (q, p) of the orthogonal system (jHj) that encounter a 
total collision at time t* and are free of diameter singularities in the 
interval [to,t*). Then there is a one-to-one correspondence between C 
and C such that, for corresponding orbits, (q, p) (t) is the orthogonal 
projection o/(q, p)(t) for every t in [to,t*). 

Proof. This equivalence follows from the fact that the orthogonal pro- 
jection is a real analytic diffeomorphism between the hemisphere z > 
of and the disk of radius k"^/^. The form of the projection proves 
the last statement of this lemma. Notice that the identification of a 
solution (q, p) of system (jlj) with a solution (q, p) of system (j9j) takes 
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place in the interval [to,t*), assumed to have no pole-geodesic con- 
figurations for (q, p) and, consequently, no diameter singularities for 
(q,p). □ 

Corollary 1. Consider a solution (q, p) of the equations of motion 
(jlj) as in the Equivalence Lemma and the corresponding solution (q, p) 
of the orthogonal system Then if c = (a,/3,7) is the constant 

vector of the total angular momentum for (q, p), the constant vector 
of the total angular momentum corresponding to (q, p) has the form 
c= (0,0,7). 

Proof. The integrals of the total angular momentum for the solution 
(q, p) of equations (jlj) are 

n n n 

{'^mi{yiZi-Ziyi),'^mi{ziXi-XiZi),^Tni{xiyi-yiXi)) = (a,/3,7). 

i=l i=l i=l 

We can write the vectors q^ and p^ of the solution (q, p) as 

qi = {Xi, yi, 0) and pi = {niiXi, rriiyi, 0), 

so the integrals of the total angular momentum for the solution (q, p) 
of the orthogonal system (Q have the form 

n 

(O, 0, cr ^ niiixiiji - UiXi)) = (0, 0, 7) =: c, 

i=l 

a fact which completes the proof. □ 



5. Total collisions 

In this section, we will use the Equivalence Lemma to generalize a 
theorem by Weierstrass and Sundman to spaces of positive constant 
curvature, first in two dimensions and then in three dimensions. 

5.1. The 2-dimensional case. We can now generalize a result known 
to Weierstrass, as communicated in a letter to Gosta Mittag-LefHer in 
1889 ([H], p. 58). Its proof was independently obtained and published 
by Sundman about two decades later, [T7|. The difference between 
this classical result and our generalization is that while in the classical 
case all three components of the constant vector of the total angular 
momentum are zero, in the curved n-body problem only one component 
must vanish. This happens because in the latter case there are no 
integrals of the center of mass, so the entire system may drift in 
before the total collapse, thus making two of the components nonzero. 
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Theorem 3. //(q, p) is a total- collision analytic solution of equations 
(jlj) as in the Equivalence Lemma, the constant total angular momentum 
vector, c = {a,P,j), has 7 = 0. 

Proof. We will use the orthogonal projection discussed in Section HI 
Since the orthogonal system describes a planar motion, the constant of 
the total angular momentum is a 3-vector with two zero components. 
The third component is exactly the constant 7 of the original system, 
as shown in Corollary [1] By the Equivalence Lemma it is enough to 
prove that, for total collision solutions of the orthogonal system, 7 = 0. 

Recall that q^ := and pi = mjCji, i = l,...,n. Then the 

orthogonal projection of the equations of motion, the energy integral, 
and the integrals of the angular momentum are obtained by replacing 
q, and pj with qj and pi, respectively. The last term of the inner prod- 
uct a ■ b = a^bx + ayby + a^b^ vanishes, so this operation becomes the 
standard 2-dimensional inner product a - b = axbx + ayby, which is well 
defined in the closed disk of radius The gradient operator is now 

the standard V with Vq^ = {dx^,dy^). Also, U^iq) remains a homoge- 
neous function of degree zero, where q = (qi, . . . , q^). Consequently, 
Euler's theorem for homogenous functions leads us to the identity (a 
proof of which can be found in ^): 

(10) q, ■ Vq,f/.(q) = 0, t = l,...,n. 

Consider now the moment of inertia /(q) = X]"=i"7,j(x^ + yf) of 
system (jlj), as defined in [1]. Notice that J(q) = /(q), where /(q) 
represents the moment of inertia after the orthogonal projection. Then 

n n 

I = 2^2 + Vl) + 2 Y.^x,dx,UM + Vidy^UM) 

i=l i=l 

n 

-2Y,miK{x^,+yl){xl + yl). 
1=1 

But from ([10]), Xidx^U^{q) + yidy^U^{q) = q^ ■ Vq^f/«,(q) = 0, so 

n 

(11) / = 2^m,(x2 + 2/^11 - <x^ + y^)]. 

i=l 

This relationship implies that / > for all t in some interval [to,t*) 
before the total collision, where to is sufficiently close to t*, therefore 
/ is increasing in this interval. We can also assume that / is either 
positive or negative, for if it changed sign at some point ti, we could 
restrict our analysis to an interval [^2, t*) with t2 > ti in which the sign 
stays the same. But if I were positive, I would increase as t — t*, so 
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the limit of I could not be zero at t*, and the total singularity would 
not take place. Therefore / < for all t in [to,t*). 
On the right hand side of the identity of Lagrange, 

mm m 2 

(12) = (j^aA) + M,-a^h)\ 

k=l k=l k=l l<k<j<m 

we ignore the first squared term, take m = 2n, a2k-i = rn]/'^Xk, Ci2k = 
mf I/., = mf i., and a,, = for = 1, . . . , n, and choose 

in the last sum only the terms that have the same index. Then the 
identity becomes the inequality 

n n n 

(13) [^mi(x,2 + y^)^ [^mi(x,2 + ^a^^j > ^^a^^^;.^. _ x.y.f_ 

i=l i=l 1=1 

But notice that 




i=l i=l 



where, obviously, p = '~f'^/n, therefore inequality f|T3|) implies that 

n 

(14) IY,m^{x',+y^)>p. 

1=1 

So if we can prove that p = 0, then •y = a ^"^j^ mi{xiyi — Xj^/j) must 
also vanish. Assume p > 0. Since, as t — t*, a total collision takes 
place, K{xf + yf) — >• for alH = 1, . . . , n, so necessarily [/^(q) — ^ oo. 
We can thus conclude from the energy relation that there exists at 
least an integer i between 1 and n for which {xj + yf) — * oo, therefore 
Yll=ii^'i + yf) ^ oo as t ^ t*. Then from identity (fTTj) . we have 

n 
i=l 

in some the interval [r, t*), with r sufficiently close to t* . This inequality 
and f|T^ imply that in the same interval we can write the inequality 

Multiplying this relationship by —21 > and integrating the ensuing 
inequality between r and t, with t in (t, t*), we obtain that 

4pln/(t) > 4pln/(r) + P{t) - P{t), 
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in which, for simphcity, we used the notation I{t) := /(q(t)). This 
relationship imphes that 

4pln/(t) > 4p\nI{T) - P{t). 

Then 

lit) > J(r)e-^'M/^'', 
an inequahty from which we can conclude that / is bounded from below 
since /(r) is finite and we assumed p > 0. But then the total collision 
cannot take place. Therefore p must be zero, and the conclusion of the 
theorem follows. □ 

Remark 3. Notice that in the proof of Theorem[3]we used only the fact 
that the potential is a homogeneous function of degree zero and that 
the equations of motion possess the integrals of the angular momentum. 
In fact the former condition is not necessary as long as one can find 
another way to show that / > in a neighborhood of t*. 

Remark 4. In the Euclidean stronger version of Theorem [3] is 

true, |14| . The proof, however, uses the integrals of the center of mass, 
which don't exist in the curved ra-body problem. 

5.2. The 3-dimensional case. We are now focusing on the curved 
n-body problem in the hemisphere 2 > of S^. With two exceptions, 
all the concepts introduced in Section [2] generalize naturally to the 3- 
dimensional case. More precisely, the position vectors of the masses rrii 
have the form = (wj, Xi, yi, Zi), i = 1, . . . , n, the gradient becomes 

and the inner product is defined as 

a ■ b := (aubu + a^K + ^yby + a^b^) 

for the vectors a = (a„, ax,ay,az) and b = {bu,bx,by,bz). Since, in Sec- 
tion [21 we expressed the kinetic energy, the force function, its gradient, 
the Hamiltonian, the equations of motion, and the integral of energy 
in terms of q and p, their formal expressions stay the same. The only 
concept we cannot naturally extend to M"^ is the cross product and, 
consequently, the integrals of the angular momentum. To bypass this 
difficulty, we will use the idea described in Section |4] of working with 
the orthogonal system instead of the original equations of motion. 

So let qj = {ui,Xi,yi) be the orthogonal projection onto the hyper- 
plane uxy of the position vector qj = (wj, x,, yi, Zi), i = 1, . . . , n. The 
vectors q^ are defined in the solid ball of radius in M"^. Equations 
(Q then describe the motion of total-collision solutions in terms of the 
orthogonal system of n bodies. The Equivalence Lemma holds in this 
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case too by replacing with S^, and the 2- dimensional disk of radius 
with the 3-dimensional ball of the same radius. 

The pole-geodesic configurations excluded from the set C of total 
collision orbits of system (jlj) occur when two or more bodies are on a 
geodesic passing through the "north pole," (0, 0, 0, k"^/^), of S^, where 
the 3-dimensional sphere is seen as a manifold embedded in M^. And 
the diameter singularities of the total-collision solutions of system iQ, 
excluded from C, show up when two or more bodies are on the same 
diameter passing through the center of the solid ball of radius k"^/^. 

Since the orthogonal system is defined in the solid ball of radius k"^/^ 
of M^, we can use the standard cross product given by 

a X b := {a,j:by - ayb^, aybu - auby, aubx - a^bu) 

for the 3-dimensional vectors a = (a„, a^,, ay) and b = 6^., by). Thus, 
we obtain for the orthogonal system Q the three integrals of the total 
angular momentum 

n 

(15) X Pi = (a,/^,7). 

To prove equations f|T5l) . notice first that 

n n n n 

Vq,t/«(q) X = ^ AijCii xqj + Y BijCii x = 0, 

1=1 i=l j = lj^i i=l 

where Aij and Bij are symmetric in i and j, the last equality following 
from the skew-symmetry of the cross product. Then, if x -multiplying 
mjCji from iQ by q^ and adding for all i from 1 to n, we obtain 

n n n 

^m^qi X qj = ^ Wn^U^^iq) x q, - ^mr^fi;(pj ■ pi)qi x q^ = 0. 

1=1 i=l i=l 

Integrals f[T^ follow by integrating the identity Yl^=i "^i^i ^ = 0. 

We can now generalize Theorem [3] to three dimensions. For the 
same reasons mentioned in Section HJ we can assume without loss of 
generality that the total collision takes place in at (0, 0, 0, 
the "north pole" of the 3-sphere. We will show that for a total collision 
solution of the original system in S^, all three constants of the total 
angular momentum of the orthogonal system must vanish. 

Theorem 4. Let (q, p) be a total- collision analytic solution of equa- 
tions (j4]) as in the Equivalence Lemma, and let (q, p) be the correspond- 
ing solution of the orthogonal system (Q . Then all three components 
a,P, and 7 of the constant total angular momentum vector belonging 
to the solution (q, p) are zero. 
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Proof. The idea of the proof is the same as for Theorem [3l Let us start 
by noticing that, according to the 3-dimensional version of the Equiv- 
alence Lemma, if (q, p) has no colhsion singularities or pole-geodesic 
configurations in the interval [tQ,t*), but experiences a total collision 
at the point (0, 0, 0, of at time t*, then (q, p) ends in a total 

collision at time t* at the origin of M^, and vice versa. We can thus 
define the moment of inertia / = ^2^=1 ''^ii'^'i + + Vi) use the 
same steps as in the proof of Theorem [3] to prove that 

n 

/ = 2 ^ 771,(^2 + + y^) [1 - /^(m^ + + y^)] , 
1=1 



in order to have the inequality / < in some interval [to, t*), with to > 
0. In the identity of Lagrange, (fT2|) . we ignore the first squared sum 
of the right hand side, take m = 3n, = fT\^'^Uk, a^k-i = T^H^Xk, 

1/2 7 1/2 . , 1/2 . , , 1/2 . f 

ask = mf! yk, 63/0-2 = Uk, hk-i = rrifj Xk, and bsk = m^! yt for 
= 1, . . . , n, and pick from the second sum of the right hand side only 
the terms that have the same index. Then the identity turns into the 
inequality 



(16) \Y^m,{uj + xj + yl 



^mi{u1 + xl + fi) 



> 



i=l 



i=l 



'^m'^^Kxiyi - XiyiY + {yiUi - ijiUif + {uiXi 



UjXi 



i=l 

But notice that 



P ■-- 



1 " 

1 = 1 



1 



i=l 



-Y^mi{uiXi - iiiXi) 

^ i=l 



iUi - yiUi) 



< 



+ 



^m,2[(xiyi - Xiyif + {yiiii - ijiUif + {uiXi - iiiXif], 



i=l 



where, obviously, p = (a^ + P"^ + 7^)/^. Thus inequality (flGll implies 
that 



/ mi{u] + x] + yf) > p. 



1=1 
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To prove that p = 0, we need to follow the same steps taken from this 
point on in the proof of Theorem [31 We are thus led to the conclusion 
that a = /5 = 7 = 0. □ 

We further show how Theorem |4] can be used to draw some conclu- 
sions about the behavior of certain solutions of the original system. 



6. Triple collisions 

In his 1912 paper, Sundman also proved that if a solution of the 
Euclidean 3-body problem ends in a triple collision, the motion must 
be planar. Again, his proof rests on the integrals of the center of mass. 

In the curved 3-body problem we might encounter solutions in which 
the triangle having the bodies at its vertices moves in without re- 
maining confined to a 2-dimensional hemisphere. But we will show 
that if the collision point is fixed in a sense we will make precise, then 
the motion must indeed take place on a 2-dimensional hemisphere. 

Theorem 5. Let (q, p) be a total- collision analytic solution of equa- 
tions (jl]) as in the Equivalence Lemma, with n = 3, and such that 
the coordinates of the corresponding solution (q, p) of the orthogonal 
system ([9]) satisfy the conditions 

3 3 3 

(17) ^niiUiit) = ^niiXiit) = ^miyi{t) = 

i=l i=l i=l 

for all t in [to, t*). Then the configuration q is confined to a 2-dimension- 
al hemisphere of curvature k > 0. 

Proof. Let (q, p) be a solution as above, i.e. one for which the coordi- 
nates of the vectors q^ = (uj, Xj, yi, Zi) satisfy the constraints uf + xj + 
yf + zf = Zi > 0, and let (q, p) be the corresponding solution of 
the orthogonal system, verifying conditions f[T71) . For system we can 
now follow Sundman's idea of proof, although the form of this system 
is different from the one of the Euclidean 3-body problem. Using (fTTjl . 
we can assume that at time t = to the three bodies of the solution (q, p) 
lie in the plane y = of M^. Since, by Theorem HJ the constant vector 
of the angular momentum is such that a = /9 = 7 = 0, we can write 
that Y^i=imiUi(to)yi(to) = YA=i^i^i{^o)yi{to) = 0, and from the last 
equation of (fTTj) . we can also conclude that X]i=i "^i2/j(^o) = 0. The al- 
gebraic system of these three equations of unknowns 2/1(^0)52/2(^0), and 
yaito) leads us to two possibilities: either ?/i(to) = 2/2(^0) = 2/3(^0) = 
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or 



det 



ni(to) 

Xiito) 

1 



1 



1 



In the former case, the solution (q, p) is confined to the plane y = 0. 
In the latter case, the bodies are initially on a straight line passing 
through the origin of the coordinate system. But we excluded the 
latter situation by eliminating diameter singularities. Thus (q, p) is a 
planar solution, therefore the position vectors of the original system are 
of the form q^ = (wj, Xj, yi, Zi), with yi = 0, Zi > 0, and uf + xf + zf = 
i = 1, 2, 3, so the bodies move on a 2-dimensional hemisphere. □ 



Remark 5. Conditions (|T71) are naturally satisfied by some classes of 
solutions of the curved n-body problem, such as the elliptic relative 
equilibria, but not satisfied by others, the hyperbolic relative equilibria 
among them. The existence of these orbits was proved in [1]. 

Remark 6. In the proof of Theorem [5] we used only the homogeneity 
of the potential, conditions (fT7|) . and the integrals of the total angular 
momentum. So this result is valid for more general potentials than the 
ones considered in this paper. 
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